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Abstract. Using the equation of motion of the damped pendulum, we in- 
troduce the averaging method on the study of periodic orbits of dynamical 
systems with small perturbation. We provide sufficient conditions for the ex- 
istence of periodic solutions of the perturbed damped pendulum with small 
oscillations having equations of motion 

9= —a0 — b0 + ef(t, 6, 9), 

where a > 0, b > and s are real parameters, with a = g/l, g the accelcration 
of the gravity, l the length of the rod and b the damped coefBcient. Here the 
parameter e is small and the smooth function / is T— periodic in t. The aver- 
aging theory provides a useful means to study dynamical systems, accessible 
to Master and PhD studcnts. 

1. Introduction 

Systems derived front the pendulum, gives to students important and practical 
examples of dynamical systems, as we can see in the case of the weight-driven 
pendulum clocks, which has had its historical and dynamical aspects studied by 
Denny in a recent paper |T| . This system has been revisited by Llibre and Teixeira 
in [2], who got the same result using averaging theory. These systems, some times, 
has also been used to introduce mathematical concepts of classical mechanics, as 
we can see in [3]. 

In this paper we will use the damped pendulum to introduce some concepts 
and techniques of averaging theory, which help us to finding periodic solutions of 
dynamical systems. 

We consider a system composed of a point mass m moving in the plane, under 
gravity foree, in which the distance between a point and m is fixed and equal to l. 
We also consider that the partiele suffers resistance to movement proportional at 
its velocity. This system is called damped pendulum. 

The position of the pendulum is determined by the angle 9 shown in Figure [1] 
The equation of motion of this system is given by 

(1) = -asin(0)-60, 

where a > 0, b > and e are real parameters, with a — g/l, g the acceleration of 
the gravity, l the length of the rod and b the damped coefficicnt. 

We consider only the motion in the vicinity of the equilibrium (9,9) — (0,0), 
i. e. we are only interested in small oscillations around this equilibrium. In this 
case, we can approximate the value of sin(0) to 9, hence, the equation ([T]) becomes 
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FlGURE 1. Pendulum. 
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FlGURE 2. Attractor Node. 



linear; moreover, denoting (x,y) = {6,6), the second order differential equation ([T]) 
linearized, can be written as the linear system of first order differential equations 



(2) 

with the following eigenvalues 



x = y, 

y = -ax - by, 



Ai = 



-b + y/b 2 - 4a 



and A? = 



-b- 



4a 



2 2 

Note that for b 2 > 4a the eigenvalues Ai and A2 are both negative, then the 
singularity (x, y) — (0, 0) is a attractor node, represented in the Figure [U Now, 
for b 2 < 4a both eigenvalues Ai and A2 have the imaginary part different of zero 
and negative real part, then the singularity (x,y) = (0,0) is a attractor focus, 
represented in the Figure [3] Both cases are topologically equivalent. If we wouldn't 
consider the damped effect, i. e. b = 0, both eigenvalues Ai and A2 would be purely 
imaginary, then the singularity (x, y) = (0, 0) would be a centre, represented in the 
Figure S) This last case is not topologically equivalent to the two cases above. 
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FlGURE 3. Attractor Focus. 
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FlGURE 4. Centre. 



To apply the averaging method, we take the damped coefficient b as a pertur- 
bation parameter of the system assuming that b = eb. Then, in this case, the 
unperturbed system is given by system ([2]) with 6 = 0, with the following eigenval- 
ues 

Ai = iuj and A2 = — iuj, 

whcre U) = y/a. 

2. STATEMENTS OF THE MAIN RESULTS 

The objective of this paper is to provide a system of equations whose simple 
zeros provide periodic solutions of the perturbed damped pendulum with equations 
of motion 

(3) x = —ax + e \—bx + f (t, x, i)) , 

where e is a small parameter and the smooth function / : R 2 — > R is T-periodic in 
t. 
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2.1. Non-autonomous perturbation. Given z = (xo,yo) G R 2 > we define the 

real functions 

(4) 

-j-. / s (bsm(2Tu) bT\ bsm 2 (Tuj) 1 f T A 
1 ^ = [ 4^ ~ T J X ° 2 J 2 2/0 ~ w / sm ( sa; )^ s ' e z ) ds ' 

_ /N Bsin 2 (Tw) / 6sin(2Tu;) 6T\ /- T , ... As . , 

J- 2 (z) = ^ ia* + f i - T )2/o + j cob(sw)/(», e As z) ds. 

A zero (xq , 2/q ) °f the system of equations 
(5) T 1 (x ,y ) = 0, J r 2 (xQ,y Q ) = 0, 

such that 

/ 9(^ 1 ,J- 2 )(z) 



dct 



z =( a; 0'2/o)/ 



9z 

is called a simple zero o f system ([5]). 

Our first main result on the periodic solutions of the damped pendulum with 
small perturbation (j3|) is the following 

Teorema 1. Assume that the function f of the damped pendulum with small per- 
turbation with eguations of motion @ is smooth and T -periodic int, then for e ^ 
sufficiently small and for every simple zero (a^jJ/o) 7^ (0)0) of the system the 
damped pendulum with small non-autonomous perturbation ([3]) has a periodic solu- 
tion x(t,e), such that (x(0, e), x(0, e)) — ¥ (ccq,j/q) when e — > 0. 

Thcorem [T] is proved in section [4] Its proof is based in the averaging theory for 
computing periodic solutions, which has been explored in section [3J 

Corolario 2. Assume that the function f of Theorem\]\has period T — (p/q)(2ir/uj), 
with p and q relatively prime positive integers. Then for e =/= sufficiently small 
and for every simple zero (xq,?/q) 7^ (0,0) of the system 

2p„ 



(G) 



Pi{z) = pTrbxo + / sin(sw)/(s, e As z) ds = 0, 
Jo 

^2(2) = Uo — / cos(su)f(s, e As z) ds = 0, 



the damped pendulum with small non-autonomous perturbation ([3]) has a periodic 
solution x(t, e), such that (x(0, e), i(0, e)) — > (xq, j/q) when e —¥ 0. 

We provide an application of Theorem [T] in the following corollary. 

Corolario 3. If f(t,x,x) = cos(u:t), then the differential equation ([3]), for e 7^ 

sufficiently small, has one periodic solution x{t, x) such that 

(x{t,e),x(t,e)) (o,l 

when e — > 0. 

The Corollaries [5] and will be proved in section [5] 
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FlGURE 5. Limit Cycle. 

2.2. Autonomous perturbation. Consider the equation of motion of the damped 
pendulum with small autonomous perturbation ([3]), with f(t,x,x) — f{x,x) given 

by 

f(x, i) — C\X + CiX + C$X 3 + C4± 3 , 

where Ci, C3 € K, C2 < b and C4 > 0. For £ ^ sufhciently small, this system has 
one periodic solution, represented in the Figure[5j 

For this system the functions !F\{x,y) and J r 2{x,y) are given by 

_ . . nx (4c 2 o; 4 — 46w 4 ) ixc\y 37rc 3 y 3 37rc 3 x 2 y 3 ■> 37rc 4 x?/ 2 
F^ X >V) = 7-^ TU 77H- + 7 1T< WX - 



Alu 5 lu 3 Alu 5 Alu 3 A Alu ' 

iry (Ac 2 lu 2 — Aku 2 ) ( 'ttc 1 i 3wc 3 y 2 \ 37rc 3 a; 3 3 2 37rc 4 t/ 3 



^(*.W) = ^ + H~ + ^ 3 ~J + ^^ + l" C4Wa;y ' 4. 

then the system of equations (JS|) has a unique solution, and this solutions is the 
trivial solution (#0,2/0) = (OjO). It implies that the Theorem[T]is not efficient in 
finding periodic solution in the autonomous perturbation case. 

Given r>0e0<^< 2ir, we define the real functions 
(7) Q(r) = —bumr + / sm(ip)f(r cos(ip), Lur sm(<p))d(p, 



n 



A zero Tq of the equation 

(8) Q(r ) = 0, 

such that dr G{rQ) 7^ 0, is called a simple zero of system ijHJ). 

Our second main result on the periodic solutions of the damped pendulum with 
small perturbation is the following 

Teorema 4. Assume that the function f of the damped pendulum with small per- 
turbation with eauations of motion ([2]) is smooth and independent of t, then for 
e =/= sufficiently small and for every simple zero rj ^ o/ the equation ([8jl, the 
damped pendulum with small autonomous perturbation @ has a periodic solution 
x(t, e), such that for all i € R, \( x (^ £ )i £ ))\ ~~ * r o> when e — > 0. 
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Thcorem 0] is proved in section |4] Its proof is based in the averaging theory for 
computing pcriodic solutions, which has been explored in section [3] 

We provide an application of Theorem [4] in the following corollary, which will be 
proved in the section [BJ 

Corolario 5. If /(z, x) = c\x + C2± + csx 3 + c^i 3 with c\ , C3 € R, C2 < b e C4 > 0, 
then the differential eguation ([3|), for e ^ sufficiently small, has one periodic 
solution x(t, x) such that 



\ {x (t,e),i(t,e))\^^^, 
when e — > 0. 

3. Averaging Theory 

We present in this section a basic result known as First Order Averaging The- 
orem. For a general introduction to averaging theory see for instance the book of 
Sanders and Verhulst [4]. 

We consider the differential equation 

(9) X = eF 1 (s,X) + e 2 R{s,X,e) 

where the smooth functions Fi : R X U — > R™ and R : R x U x (—e f, s f) — > R™ are 
T-periodic in the first variable and U is an open subsef of R". 
We define the averaged system associated with system © as 

(10) Y(t)=Eh(Y(s)), 
where /1 : U — > R™ is given by 

(11) h{Z)= [ F 1 (s,Z)ds. 

Jo 

We present, in the following theorem, the conditions for the singularities, of the 
system (fT0|) , be associated with the periodic solutions of the differential system ((9j . 

Teorema 6. Assume that for a G U with /i(a) = 0, there exist a neighborhood V 
ofa such that f x (Z) ^ for all Z e V - {a} and det(<#i(a)) ^ 0. Then, for \e\ > 
sufficiently small, there exist a T-periodic solution X(t,e) of the system ([9]) such 
that X(0, e) — > a as e — > 0. 

Theorem |5] can be announced with weaker hypothesis using Brower degree. For 
a proof of Theorem [S] see [S] . 

4. Proof of Theorem [1] 

In order to use the Theorem [6] in the proof of Theorem [TJ we have to modify the 
equation ([3]). Let's take y = x. If we denote 

y)' A= (- 0) - d F ^={~- b y + °f(t, 
then the equation ([3]) can be written in a matrix form 

(12) ± = Ax + eF(t,x). 
We define a new variable y (t) € R 2 by 

y(t) - e- Ai x(i), 
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where e At is the matrix of fundamental solution of the unperturbed differential 
system (p2|), Le., e = 0. Note that y(0) = x(0). 

The system (|12p is written in the new variable y as 



(13) y = ee- At F(t,y), 

with F(t,y) = F(t,e At y)- 

Now, we are ready to prove the Theorem [1] 

Proof of Theorem\7\ Assume that the smooth function /(f, x) is T-periodic in t, 
then the smooth function F(t, y) is also T-periodic in t. We shall apply Theorem 
|6]of the section[3]to the differential equation (fT3|). We note that equation (fT3| can 
be written as equation taking 

X = y, a = t, F l {t,y)=e- At F(t,y) and R{t,y,e) = 0. 
Given z = (:co, yo) G R 2 , we define the averaged function fi : M 2 — > R 2 

A (z) = / F x (s,z)d8. 
Jo 



Proceeding with calculations we have 



, . fesin(2Tw) bT\ bsin 2 (Tuj) f 1 sin(sw)/(s, e As z) , 



6sin 2 (Ta;) / 5sin(2Tu;) bT\ f T , , 

— 2 — + l — ~ T ' yo + ' cos ( sw )/( s ' e z ) ds 

(Ji^TMz)). 



n 



Then, by Theorem |6] we have that for every simple zero (xq,j/o) S V' of the 
system ([5]), there exist a periodic solution y (t, e) of the system (fT5)) such that 

y(0,e) (x5)2/o) when e ^ 0. 

It has proved the existence of a periodic solution x(£, e) = (x(i, s), y (t, e)) of the 
system (fT2|) and consequently of the equation (j3|). As x(0) = y(0) it follows that 

x(0,e)\ (x* Q 

y(o,e)J U* 

when e — > 0. Hence Theorem [T]is proved. 

□ 

5. Proof of Theorem [4] 
Rescaling the time by r = uit, the equation becomes 
(14) x" = -x + e f (t, x, x') 

where now the prime denotes derivative with respect to the new time t and 

/ (t, x,x') = — (—bojx' + f (t, x, ux')) . 

a v ' 
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In order to use the Theorem [6] in the proof of Theorem |4j we have to modify the 
equation (HJ). Take y = x' , then the second order diffcrential equation (fT4|) can be 
written as the system os first order differential equations 



(15) 

In polar coordinates 



z' = y, 



y' = -x + ef(t,x,y). 



x = r cos(yj), r 2 — x 2 + y 2 , 

y 

y = rsin(<p), t<m((p) = -, 

x 

with r>0e0<(/3<27r, the system (fT5"]) becomes 

r' = eaw((p)g(r,<p), 

(16) 

(p' = -1 + e y ' x ' , 

r 

where g(r,(p) — f (r cos((p) , r sm(tp)) . 

Using ip as the new independent variable, the system (fTBf becomes 

dr er sm(tp)g(r, ip) 
dtp e cos(ip)g(r, cp) — r 

Expanding the above expression in Taylor's series near e = 0, we have 



(17) — = -e 8w{tp)g{r, ip) - e tan(p) ^ e 

' n— 2 



2 



Now, we are ready to prove the Theorem 2] 

Proof of Theorem^ Assume that the smooth function f (t, x, i) doesn't depend of 
the variable t, i. e, f (t, x, x) = f(x, i), then the smooth function g(r, ip) also doesn't 
depend of the variable t t. We shall apply Theorem [6] of the section [3l We note 
that the equation (jTTJ) can be written as equation © taking 

X = r, s = (p, Fi(ip, r) = — sxD.(np)g(r, (p) and 



ra-2 



Observe that the function F\((p,r) is 27r-periodic in ip, then the hypothesis, about 
equation ([9]) of Theorem [6j are satisfied. 



Given r > 0, we define the averaged function / : 



/•27T 

h(r)= / F x (<p,r)dtp. 
Jo 
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Proceeding with calculations we have 

/-27T 

/i(r) = / am(<p)g(r, ip)dip, 



2tt 

sin(p) (/(r cos(<£>), r sin(p))) dip, 



o 

2-7T 



sin(y)— (— !wr sin(yj) + /(r cos(^j), wr sin(y))) d^, 
sin(<p)(r cos(</?), wr sin(</?))d</j, 



buj r .2/ u 1 r 

— r / sm ((p)d(p / 

a 7o a J o 

1 / sin<p(r cos(ip),u)rsm(tp))d(p, 



a a _ 

Thcn, by Thcorcm|6l we have that for every simple zero Tq for the equation (|8|). 
there exist a periodic solution r(ip,e) of the equation (fl7|) such that 

r(0, e) — > Tq quando e — > 0. 

It has proved the existence of a periodic solution jc(t, e) of the equation ([3]) such 
that 

|(i(i,e),a;(f,e))| rS, 
for all t G R, when e — >• 0. Hence Theorem|3]is proved. □ 

6. Proofs of Corollaries 

Demonstragao do Coroldrio\^ We have that f (t, x, x) is T-periodic with T — (p/q)(2n/u), 
then T = qT is also a period of the function f (t, x, x), Since sin(2Tw) = 0, it follows 

i-T 



Note that the systems © is equivalent to the systcm 

Ji(z) = 0, J" 2 (z)=0, 

because both equations only differs in a non-zero multiplicative constant. Hence, 
applying the Theorem[TJ the result of the corollary [2] follows. □ 

Proof of Corollary^ We have that f(t,x,x) is T-periodic, with T — 2tt/lu, then 
we can apply the Corollary [2j Proceeding with the calculations we get 

■£■ r \ bn 
-O (z) = — x , 

U! 

~ bn 7T 

Tiyz.) = —y — , 

U! UJ 

with 

(«) = k = 



bT 




Y X0 ~ 


-S 




u J o 


bT 






f' 
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as a solution of the system (j6]). Sincc 



det 



dz, I ~ 



2 



it follows that (xq, j/q) ls a simple zero of system ([6]). Hence, applying the Corollary 
[2j the result of the corollary [3] follows. □ 

Proof of Corollary [31 Proceeding with calculations we get 

3 

G (r) = —oj 3 nc4r 3 — (b — C2)uinr, 
a polynomial in r, whose zeros are 



and ± 2 lb - C2 



uj y 3c4 
Let's take 



rS=2 J'- 02 



V 0C4 
since 

j-{r$) = 2wtt(&- 02)^0, 

it follows that r g is a simple zero of system ([8]l, Hence, applying the Theorem[4l 
the result of the corollary [5] follows. □ 



7. CONCLUSIONS 

The averaging theory is a collection of techniques to study the persistence of 
periodic orbits of a dynamical system under small perturbation. 

In this paper, we present some of this techniques and we use it to give conditions 
for the existence of periodic solutions of the equation of motion of the perturbed 
dampcd pendulum. This problem is divided in two cases: non-autonomous pertur- 
bation and autonomous perturbation. The same technique is used in both cases, 
diffcring only in the preliminary treatment in each case before apply the first order 
averaging theorem. 

For the first case, the equation of motion is written as a system of two first order 
differential equation, which the independent variable is the time. For the second 
case, using polar coordinates, we approximate the equation of motion by a first 
order differential equation independent of the time. 

We get, in both cases, systems of equations that have its simple zeros associ- 
ated with periodic solutions of the equation of motion of the perturbed damped 
pendulum. These relations are explicit in the theorems that we prove here. 
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